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Abstract. We here specialize the standard matrix- valued polynomial inter- 
polation to the case where on the imaginary axis the interpolating polyno- 
mials admit various symmetries: Positive semidefinite, Skew-Hermitian, J- 
Hermitian, Hamiltonian and others. 

The procedure is comprized of three stages, illustrated through the case 
where on iR the interpolating polynomials are to be positive semidefinite. 
We first, on the expense of doubling the degree, obtain a minimal degree 
interpolating polynomial P(s) which on iR is Hermitian. Then we find all 
polynomials ^(s), vanishing at the interpolation points which are positive 
semidefinite on iR. Finally, using the fact that the set of positive semidefinite 
matrices is a convex subcone of Hermitian matrices, one can compute the 
minimal scalar /3 > so that P(s)+/3^!(s) satisfies all interpolation constraints 
for all > 0. 

This approach is then adapted to cases when the family of interpolating 
polynomials is not convex. Whenever convex, we parameterize all minimal 
degree interpolating polynomials. 



1. INTRODUCTION 

Probably the simplest version of interpolation problem is as follows. Given a family 
of functions J-, nodes x\ , ... , x p and image points Y\ , ... , Y p , search for F G T 
so that 

(1.1) Yj = F( Xj ) .7 = 1 , ... , p. 

More specifically, find out whether such F(s) exists and if yes, search for all "simple" 
interpolating functions in J- ' . In the context of rational functions, "simple" means 
low degree, which here takes the form of the McMillan degree, see e.g. [15], [18] . 

There is a vast literature on this classical problem. For a comprehensive study see 
e.g. [18]. Additional relevant literature is presented in Section [2] 

Using the framework in (|1.1[) . in this work, the nodes Xj are in C, the image points 
Yj are m x m matrices and T are polynomials F of a complex variable s, i.e. 

9 

(1.2) F(s) = J2 °ks k C fc e C mxm . 

k=0 



1991 Mathematics Subject Classification. 11C99, 15A99, 32E30, 41A05, 47A57, 47B65. 

Key words and phrases, interpolation, polynomial, matrix-valued, structured matrices, La- 
grange, generalized positive functions, generalized positive even, convex cones, convex invertible 
cones. 

This research was supported in part by the Bi-national Science Foundation grant 2010117. 
D. Alpay thanks the Earl Katz family for endowing the chair which supported his research. 

1 



2 



D. ALPAY AND I. LEWKOWICZ 



Recall, that the McMillan degree of matrix- valued polynomials is well defined. In 
[20l Corollary 2.1.1] it is shown that for F(s) in (|1.2|) it is equal to the rank of the 
block-triangular, block- Toeplitz matrix 
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In particular, if C q is nonsingular, then the McMillan degree of F(s) is equal to 
m(q + 1). 

In this work, the polynomials in J- are restricted to have various symmetries on 
the imaginary axis, described in the sequel. The importance of polynomial matrix 
interpolation with symmetries, was raised in [331 Subsection 2.1.48]. In the frame- 
work of real variable, this problem has already been treated in |35) . That work 
differs from ours in many ways. 

Through Examples ll.il HOI 12. 2[ 16.31 and 17. ll part A, we illustrate the fact that even 
in the scalar case, the question addressed here is not trivial. We start with the 
following. 

Example 1.1. As a prototype example we shall seek a minimal degree interpolating 
function F(s) so that 

F(s) (2) -> (ft 



, 3 J V 50 / 

Here and in Examplcs l2.1[|2~2"l and l6.3l we shall consider various families of functions 
J 7 with the same points. 

Taking the family J- to be unstructured polynomials a straightforward computation 
yields, 

Fi(s) = -45s 2 + 180s- 121. 

□ 

As mentioned, we here focus on structured polynomials. To formally set-up the 
problem addressed here, we need some background. 

1.1. Functions with symmetry on the imaginary axis. Let F be m x m- valued 
rational function in the sense that 

F:C -> C mxm . 

In the sequel we shall use the following notation, 

F*(s) :=F*(-s*). 

We shall denote by C+ (C+) the open (closed) right half plane and by F m (P m ) 
the sets of m x to positive semidcfinitc (definite) matrices. Whenever clear from 
the context, the subscript will be omitted and we shall simply write P or P. 

We call functions F(s) Positive, denoted by V, if they are analytic in C+ and 

(1.4) (F{s) + F*{s)) eP se<C+. 

These functions has played an important in the theory of electrical networks from 
around 1930, see e.g. [15] and [22] . They also serve as the corner stone of the 
theory of linear dissipative systems (a.k.a absolutely stable), see e.g. [HI Theorem 
2.7.1], [221 3.18]. 
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One can relax the condition and call functions F(s) Generalized Positive, denoted 
by QV, if 

(F(s) + F*(s)) e P almost for all s e iK. 

For an early study of rational QV functions see [2] and for recent references see 
e.g. [5], [IH] and [TT]. These functions were studied in different frameworks, for 
example when the upper half plane replaces C+, see e.g. [33], [3(5], [37] or when the 
unit disk replaces C+, either in the domain see e.g. [30] or the image, see e.g. [39] . 

Abusing terminology of real scalar functions, we call a function F(s) Odd if 

(1.5) F*(s) = -F(s). 

This implies that on the imaginary axis F(s) is skew-Hermitian, i.e. 

We shall denote by Odd the set of odd functions. Note that Odd C QV, for details, 
see [TT] Proposition 4.2]. 

In a similar way we call F(s) Even if 

(1.6) F*(s) = F(s). 

This implies that on the imaginary axis F(s) is Hermitian, i.e. 

(F(sh^)* = F(s) Um . 
The set of even functions, denote by Even, was studied in |11[ Section 5], 
Of particular interest is the class of Generalized Positive Even, denoted by QVE, 

/j 7) F#( s )=F(s) 

F(s)£F almost for all s£iR. 

For details see [TTJ Section 5]. Recall that F € QVE if and only if there exist G(s) 
so that 

(1.8) F(s)=G(s)G#(s). 

If F(s) is analytic on the imaginary axis (|1.8p is called spectral factorizatioi^\, see 
e.g. PI] Section 5.2], [UJ Chapter 9], and [H Section 19.3] Else, (JTSJI is a pseudo 
spectral factorization, see e.g. |21i Chapter 10], 

For future reference we recall that a convex cone which in addition is closed under 
inversion is called a Convex Invertible Cone, cic in shor1@ , see e.g. [315], [37J and 

It is easy to see that the sets Odd and Even are closed under positive scaling, 
summation and inversion, i.e. cics. 

The set V is a subcic of QV functions. More precisely, V is a maximal cic of 
functions which are analytic in C+, see e.g. [55] Proposition 4.1.1]. The cic 
structure of QV functions was studied in (TTJ . 

Recall that the intersection of cics is a cic, e.g. [26] Proposition 2.2]. In particular 
QV£ is the intersection between the cics of QV and Even. We now describe an 



^G(s) is analytic in C+ sometimes also G~ 1 (s) is analytic there. 

2 Strictly speaking, this means that whenever the inverse exists, it also belongs to the set, e.g. 
the set of positive semidefinite matrices is a cic. In contrast, the open upper half of C is not. 



4 



D. ALPAY AND I. LEWKOWICZ 



intermediate set between the Even its subcic QP£. Recall that Even, see (|1.6[) . 
means F\ m is Hermitian and QP£ means F| ; is positive semidefinite, see fj 1 . Tj) . 

Next consider the set where for all 10 G K, F(iuj) is Hermitian and there is no 
eigenvalues crossing from C_ to C+ (or vice versa). Roughly, along the imaginary 
axis the inertia is (almost) fixed. More precisely, we call a m x m-valued ratio- 
nal function F(s) v -Generalized Positive Even, denoted by vQP£ if it admits a 
factorization 

F(s) = G( S )diag{-/„ , I m ^}G#(s) v G [0, l], 

for some G(s). In particular, for v = one returns to the QP£ case in (|1.7|l and if 
v = r7i, then — F G QP£. Thus, to avoid triviality, in the sequel we shall focus on 
vQP£ functions admitting factorization 

(1.9) F(s) = G(s)diag{-/„ , 7 m _„}G#(s) v G [1, m - 1]. 

If F(s) is analytic on the imaginary axis (|1.9|) is called J-spectral factorization. 
Else, this is ./-pseudo spectral factorization, see e.g. [21] Part VII]. 

It should be pointed out that the technique wc use in the sequel does not require 
factorization of vQP£ (or QP£) functions. 

The aim of this work is to offer a way for solving the interpolation problem in 
where the family F is comprized of matrix-valued polynomials within: Odd, Even, 
QVZ andvQP£. Moreover, whenever F is convex, all minimal degree interpolating 
polynomials are given. 

We now set up the main idea of this work in terms of (jl.ip , through the framework 
of F = QP£ polynomials. Other cases will turn to be (not necessarily small) 
variations on the same theme. 

(i) In classical unstructured polynomial interpolation one obtains P(s) (P : C — > C mxm ) 

q-l 

(1.10) P(s) = J2 C kS k C fe GC mxm , 

k=0 

an unstructured interpolating polynomial of minimal degree, namely it is 
at most mq with p — 1 > q (recall p is the number of interpolating nodes) . 
On the expense of doubling its degree (i.e. 2m(77 — 1) > deg(P)), one can 
introduce partial structure to the interpolating polynomial P(s) in (jl.lOp . 
see Proposition 13. 11 In particular, the resulting polynomial P(s) may be in 
Even, see Subsection 15.11 or in Odd, see Subsection 15.31 
If the obtained P(s) is already in QPS , we are done. Assume that this is 
not the case. 

(ii) It is easy to obtain all minimal degree QPS polynomials vanishing at 
x\ , ... , x n , denoted by ^(s), (2?7m > deg(^)), see Section [3] and Propo- 
sition [BTTJ 

(iii) Using the above P(s) and 'J'(s), for all f3 G C, 

(1.11) F(s)=P(s) + f3^(s), 

is an interpolating polynomial. In fact, F(s) is in Even for all /3 G M. 
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(iv) The family of QV£ polynomials is a convex subcone of Even and by con- 
struction deg(^E') > deg(P). Thus, one can find f3 so that for all /3 > $ the 
interpolating polynomial F(s) in ([1.11)1 is in QVE . See Proposition ^. 21 

As mentioned, interpolation with QVE polynomials serves as a prototype of our 
technique. In Section [5] the same idea is extended to polynomials with various 
symmetries on the imaginary axis. In Section [7] refine the procedure by exploiting 
possible special structure of P(s) obtained in stage (i) of the "recipe" to enlarge 
the family of polynomials ^f(s) obtained in stage (ii) of the "recipe" and to reduce 
the minimal degree of the interpolating polynomial F(s). In Section [8] we modify 
the recipe to allow interpolation with non-convex set of polynomials, e.g. vQVE in 
(|1.9p . Finally, in Section [9] a sample of future research problems. 



2. MOTIVATION AND BACKGROUND 

In this section we review various aspects of the interpolation problem at hand. For 
completeness we start by presenting three of the popular variants of the interpola- 
tion problem not addressed in this work. 

First, the tangential interpolation, i.e. given x\, ... , i p S C, Hi, ... , v p G C m 
and Yi, ... , Y p G C l find F(s) in T (F : C -> C' xm ) so that 

(2.1) Y s = Ffafa j = l, ... , p. 

Strictly speaking, this is the right tangential interpolation problem. The left tan- 
gential interpolation problem is: For given x\, ... , x p G C, u\, ... , u n G C l and 
Z t , ... , Z p G C m finding F(s) (F : C ->■ C' xm ) so that 

(2.2) Z; = u*F( Xj ) j = l, ... , p. 

Combining both, one obtains the bi-tangential interpolation problem. For sample 
references see e.g. [20], [38] for the polynomial case, [16] for the unstructured 
rational case, and [4], [5], for the structured rational case. 

As a second popular variant we mention that in control theory there has been 
an interest in problems of the following form. Given A(s),B(s) m x m-valued 
polynomials where in addition B G Even one seeks all m x m-valued polynomials 
X(s) satisfying, 

X*(s)A(s) + A*(s)X{s) = 2B{s). 
See e.g. [41] and the survey in [31] 

To present a third popular variant, we start by resorting to the notion of re- 
verse polynomial. Recall that F(s) is said to be the reverse of a polynomial 

F(s) = J2 s J C 3 in if 

(2.3) F{s):=s^F(s- l ) = J2s j C q - j , 

3=0 

see e.g. [H Section 2], [33 Eq. (8.17)]. 

Next, we recall in the notion of linearization. For constant mq x mq matrices A, B 
we say that the pencil sA — B is a linearlization of the polynomial F(s) in (|1.2p if 
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there exist mq x mq- valued unimodular polynomials E(s), -D(s), G(s), H(s) so that, 

diag{F(s), I m{q - X) } = E{s){sA-B)D{s) 
diag{^( S ), /„,(,_!)} = G(s)(A-sB)H(s). 

For details see e.g. [15] . 

Roughly speaking, out of the many variants of interpolations, two of the better 
studied frameworks are (i) unstructured polynomials of the Lagrange type see e.g. 
[13] . [23] . [45] . [48] and (ii) structured rational functions, of the Nevanlinna-Pick 
type, see e.g. [T5] Section 18]. 

There are fundamental differences between these two problems: The Nevanlinna- 
Pick type problems is far more involved, but the underlying structure allows for 
the use of powerful tools. 

Addressing interpolation through structured polynomials, as we do here, turn most 
of the classical arsenal, like linear fractional transformation, irrelevant. 

Nevanlinna-Pick interpolation of positive functions V has been well studied, see e.g. 
[TBI Chapter 18]. It was extended, not in the framework of QV but of (i) generalized 
Schur functions (contractive on the unit circle), in numerous works, see e.g. pQ, 
[T7] , [H] , [25] , and [53] and (ii) generalized Nevanlinna functions (mapping the real 
axis to the upper half plane) [2J, [12], [7J, [M], E3 Section 3] and [BJ. 

It should be pointed out that this extension of Nevanlinna-Pick interpolation from 
addressing V to QV is computationally involved. Moreover, the existing parameter- 
ization of all QV interpolating functions (after being translated from the generalized 
Schur framework) neither single out QVZ functions nor polynomials. This is illus- 
trated next. 

Example 2.1. Consider the interpolation points in Example ll.il 

I. Assume now that the family T is the set of rational QV functions. Then, minimal 
degree interpolating function is 

150(574 - 451s) 
2 ^ ~ 41(66 -41s) ' 

Clearly, Fi{s) is neither in Even nor a polynomial. 

II. Consider the interpolating polynomial in Example 11.11 Then, F\{s) is neither 
in Even nor in QV. □ 

To summarize, from a practical point of view, the existing interpolation scheme for 
QV functions, is not very helpful for interpolation by QVE polynomials. 

The Lagrange approach 

A classical approach, see e.g. [57] Section 2.10], common to some interpolation 
problems of the form (|1.1[) is here illustrated 

Take, 

v 

F(s) = Y,*j(*) J", j = l, ... , p 

3=1 
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so that 

XI X2 ... X p 

#i(a) Y 1 ... 
F 2 (s) Y-2 ... 

F p (s) ... Y p . 

In the framework of scalar (non-structured) polynomial interpolation, this approach 
probably preceded Lagrange [55] (in [IS] it is attributed to |50j). 

Computationally motivated, an interesting choice of Fj(s) for unstructured scalar 
polynomials was presented in [331 Eq. (4.2)]. 

The above straightforward approach to interpolation problems has some limitations: 

(i) It first assumes that it is easy to construct the elements Fj in J- (for example if 
T is the set of QV polynomials, this is not easy, see e.g. [TT| Example 5.3b]). 

(ii) It assumes that J 7 , the family of interpolating functions, is convex (in Section 
|SJ we address a non-convex family J 7 ) . 

(iii) In addition if for example J 7 is a convex set of rational functions this scheme 
may yield high degree interpolating functions. 

We now further scrutinize this interpolation scheme. To this end we here employ 
it to an exampl^f). 

Example 2.2. We here show that adapting the Lagrange approach to the case 
where J- = QV£ polynomials, enables us to construct some of the minimal de- 
gree interpolating polynomials in the problem addressed in Examples 11.11 and 12.11 
Indeed, take 

F 3 (s) = Fi(s) + F 2 (s) + F 3 (s) 

with 

#!(«) = (4- S 2 )(9-s 2 )(Q(l-s 2 ) + f) a>0, 
F a (a) = (l-s 2 )(9-s 2 )(/3(4- ;i 2 )-| ;s 2 ) ^>0, 
F 3 (s) = (l- s 2 )(4-s 2 )(7(9- s 2 ) + f ) 7>0. 

This can be aggregated to a single parameter 

(2.4) F 3 (s) = -3s 4 + 34s 2 -13 + /?(l-s 2 )(4-s 2 )(9-s 2 ) 13 -\ ■ 

□ 

In Example 16.31 below we show that in contrast to to (|2.4I) . there is a QV£ interpo- 
lating polynomial for all > 4 . Thus, adapting the Lagrange approach to QV£ 
polynomials is appealing due to its simplicity. However, even in the scalar case it 
turns out to provide conservative results. See also Part A of Example 17.11 

We conclude this section by further examining the structure of some the families J 7 
involved. Let F(s) be a m x m- valued polynomial as in (jl.2p . It is straightforward 
to verify that 

J~ C'2k C2fc+1 

Even Hermitian skew — Hermitian 

Odd skew— Hermitian Hermitian. 

However that there is no explicit way to characterize the sets QV, QV£ vQV£ only 
through the structure of their coefficients. For example, scalar, second degree, QV£ 



^We have already used it in the framework of QVS polynomials in 1111 Example 5.3a]. 
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polynomials are given by 

a + b(r + is) 2 with a > 6 > rel 

(it is a strict subset of all F(s) in (jl.2p with q = 2 where C > 0, Ci € iK and 
> C2). Hence, it is conceivable to presume that to solve the problem in (jl.ip one 
needs to go beyond a simple modification of the classical polynomial interpolation. 



3. PARTLY STRUCTURED POLYNOMIAL INTERPOLATION 

The classical to x m-valued unstructured polynomial interpolation can be formu- 
lated as follows. Given x±, ... , i p eC (distinct) and Y\ , . . . , Y p £ C mxm find 
P(s), a minimal degree polynomial (P : C — > C mxm ) see (jl.lOp , so that 

Y j =P{x j ) j = l, ... , p. 

It is known that the problem is solvable and in (jl.lOp , p > q. 

We here adaptQ an idea from [4], [5j Section 2] enabling us, by roughly doubling q 
in (jl.lOp . to impose structure on the matrical coefficients C'k- The problem can be 
formulated as follows. 

Given: xi, ... , x p G C and A, B, F x , ... , Y p <= C mxm find P(s) a low degree 
m x to- valued polynomial (P : C — > C mxm ) so that 

. . P*(s) - AP(s)B A,BeC mxm non-singular 

P(x 3 ) = Y } j=l, ... , p. 

To guarantee feasibility of the problem one needs to assume that the data satisfies 

,„ „. xj-x k =0 =>• Yj = Y k 

(3.2) , p > k > j > 1. 

Next, we shall call the interpolation data set reduced if out of feasible points 
X\ , ... , x p we extract a maximal subset x\ , ... , x n , i.e. p > n (and the 
corresponding Y\ , ... , Y n ) so that 

(3.3) I n > k > j > 1. 

Note that the choice of x\ , . . . , x n out of the p given nodes, is not unique (n is 
unique) but it will not affect the proposed procedure belowQ. 

Proposition 3.1. Given a feasible data set (|3.2p . There always exists an interpo- 
lating polynomial P{s) (jl.lOp satisfying (|3.ip with 

2n > q, 

where n is the dimension of the reduced data set (|3.3p . 

The coefficients C , ... , G-in—X ar £ obtained from the following matrix equation 

(3.4) XC= Y 



^Originally, it appeared in other frameworks: Caratheodory functions in [4] and the unit disk 
H2 functions in [5]. 

^To summarize, in X mis the dimension of F{s), the number of given interpolation points is p 
and it is then reduced to n. Finally q is the number of matrical coefficients in Ql.100 and 111, lift . 
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where the dimensions of both C and Y is 2nm x m, 

( Yl \ 



C:= 



C„ 



,c 2 „ 



Y:= 



(AY 1 B)' 



\ {AY n B)' I 

and X is the 2nm x 2nm block- Vandermonde matrix 



x= 



Xilfj 



Im X„I m </ m 
Im ~X*I m (-xl) 2 I„ 
I m -x*I m (-xi,) 2 I„ 



\I m -x" n l m (-x* n ) 2 I m ... (-x n *) 2n - 1 I m I 

Proof Note that X can be written as 

X = X®I m , 

where X is an actual 2n x 2n Vandermonde matrix (i.e. substitute m = 1 in X 
and (g> denotes the Kronecker product, see e.g. [HI Section 4.2]. Minimality of the 



data implies that the 2n points x\, 



-x* are distinct, so X is 



non-singular, see e.g. [44l Exercise 12]. This in turn implies that X is nonsingular, 



see e.g. [21 Corollary 4.2.11]. Hence, the matrical coefficients C , 
([l.lOj) are unique and explicitly obtained. 



C 2n -i in 



k = 0, 



Finally, the structure guarantees that 
C k * = {-l) k AC k B 



In the sequel we focus on a special case of (|3.ip where 

B = (A*)- 1 

Namely P(s) is so that, 



2n- 1 



□ 



(3.5) 



P#(s) = AP(s)(A*)- 
P{x 3 ) = Y 3 



non- singular 



j=l, ... , n. 



4. CONSTRUCTING MINIMAL DEGREE SYMMETRIC NEUTRAL POLYNOMIALS 

We now address the problem of constructing, within a prescribed family J 7 , poly- 
nomials $(s) of minimal degree, vanishing at the given nodes x\, ... , i p £C. 

Using (|3.3p let now x±, ... , i n be a resulting reduced set, see (|3.3I) . Next, take 

n 

(4.1) *(s) := Y[( Xj -s)M( X ;-+s) MeC mxm , 

i=i 

with M parameter. Clearly, ^(s) vanishes at the original points x\, ... , x p . 

Hence, taking M to be: Hermitian or Skew-Hermitian yields a minimal degree \1/ 
in: Even or Odd, respectively. 
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Note that for s = iuj with ui real, *(s) in (|4.1[) satisfies 



(4.2) V(iu)) = M J| 



Xj — IUJ 



2 



To guarantee that the sets of polynomials described in (|3.5[) and (|4.1I) indeed in- 
tersect, we need to further restrict the parameter M in (|4.1[) so that the product 
AM, with A from Q3.5p . is Hermitian, i.e. 



(4.3) AM = (AM)* 



5. INTERESTING SPECIAL CASES 

We here specialize the "recipe" from Section [1] to interesting classes of polynomials. 

5.1. Even. To guarantee feasibility of the problem, one needs to substitute in fl3.2[) 
A = B = I m and thus obtain, 

xj-x k =0 Yj = Yfc 

p > k > j > 1. 

Taking in (|3.5[) A = J m results in P(s) is in Sven, i.e. P = 

5.2. J — Sven. We shall find it convenient to denote by J an arbitrary m x m 
Hermitian involution, i.e. 

(5.1) J = J* = J-\ 

Recall that J is unitarily similar to diag{— , I m -v} with v 6 [0, m], see e.g. [43l 
Theorem 4.1.5]. 

Taking in (|3.5[) A = J with J as in (|5.1[) . results in P(s) in J — Even. 

In order to have the problem feasible one needs to assume that the original data 
satisfies 

Xj+xl=Q =s. y 3 ■ = jy'j p>k>j>l. 

i o -I m \ 

For example if m is even and J = then on the imaginary axis P(s) 

W y 

has Hamiltonian structure, see e.g. (35J Section I], [40l Problem 3.21, Eq. (6.3.3), 
Theorem 11.5.1] or for a thorough treatment of this structure, [46l Section 7.2]. 

5.3. Odd. Taking in p.5[) A = il m results in P(s) in Odd. In particular, in order 
to have the problem feasible one needs to assume that the original data satisfies 

Xj +x* k =o =^ Yj = -y; P>k>j> 1. 

In this case, the condition in (|4.3|) implies that M in (|4.1[) is skew-Hermitian. 
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6. MINIMAL DEGREE INTERPOLATING GENERALIZED POSITIVE EVEN 

POLYNOMIALS 

Let x\ , ... , x p <G C and Yi, . . . , Y p <E C mxm be a feasible data set i.e. 

Xj+xl=Q Y j= Y k p>k>j>l 

Xj GiR =>■ Y,GP m j=l, ... , p. 

One searches all minimal degree interpolating QP£ polynomials, F(s) in To 
simplify presentation assume that the data is already reduced, i.e. 

p = n. 

Now take: (i) From Subsection 15.11 a minimal degree interpolating P € Even, (ii) 
From (|4.ip all minimal degree "J G 5'P£ vanishing at the interpolation points: 

n 

(6.1) *(«) := J]^ - s)M(x; + s) M £ P m . 

We now establish the minimality of the degree of *ff(s) in (|6.1|) . 

Proposition 6.1. Let \l/(s) be a m x m-valued, minimal degree QV£ polynomial 
vanishing at a given set of distinct points x\, ... , x p € C. Le< a/so xi, ... , i„6C 
&e a corresponding reduced set, see (|3 . 3[) . XTien vE'(s) is o/ the form (|6.ip . 

Proof : From a given data points Xx, ... , x p let us denote 

p 

g(s) ■= -s). 

Recall that every QV£ function admits a factorization of the form (|1.8|) . Thus, 
every m x m-valued C?'P£ polynomial ^(s) vanishing at of the form 

y(s) = g(s)y(s)g*(s) ^ e GPS. 

Let now x± , ... , x n be a reduced subset of the data points and then 

n 

9o(s) := ]J(xj - s). 

i=i 

Hence, without loss of generality one can take 

9(8) = g (8)*(8)gf(a) VeGPS. 
To guarantee minimality of the degree of \P(s), take it to be of degree zero, i.e. 

*(«) = M g P m , 

so the claim is established. □ 

Without loss of generality we shall find it convenient to normalize M in (|6.1I) so 
that 

l|M|| 2 = l. 

Following the "recipe" from Section [TJ the sought interpolating polynomial F(s) is 
of the form 

(6.2) F(s)=P(s)+/3*(s), 

with P(s), ^(s) from Propositions 13. ll!6. 11 respectively and /3 a parameter. 
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Assume that P(s) in (fTTTOl) is not in QVE. First, we fix in (|6.1j) an arbitrary M in 
P m (i.e. non-singular). By construction, for all (3 £ R in ()6.2[) is an interpolating 
polynomial (|l.ip in Even. 

Now, on the one hand in (1 1 . 10[) P £ Even is of degree of (at most) (2n — l)m. 
On the other hand, in (|6.1I) ^f(s) is of degree 2nm and in QV£ , a convex subset of 
£uen. Thus, there exists j3 > so that in f|l . 1 If) . ([Oil 



* e g-p£ V/3 > p. 

We next show that here $ can be explicitly obtained. 



Proposition 6.2. Assume that P(s) in (| 1 . 10[) is not in QVE and that in (|6.1 
M £ P m is given. Then j3 in (jl.llj) . (|6.2|) is given by 



3.3) /3 



mm 



n 



xj — ilu\ mm 

i— 1 , ... . 



A, Af" 1 ^ i k C k uj k 



k=0 



Proof Indeed, note that F(s) in (| 1 . 1 1 1) . (|6.2|) is in QVE if the following relations 
hold 



for almost all s£il 
for almost all cj^IE 

for almost all cjGI 



F( S ) 
F(i«) 

fc = j = i 



6F„ 



'2n — 1 1 1 n 

for almost all wgl I £ M ~ 5 c fc M ~ 5 icj fc +/3 f] ) 
fc=o 3 = 1 



for almost all w£lE 

This can be written as, 



2n— 1 _1 _I 
E M 2 Cj.M 2 , 
fc=0 

n i^-m 2 

3 = 1 



6P„ 



for almost all LJEff 



for almost all w( 



for almost all ufEff 



for almost all a;£lR| FJ | £c j — icj 

Thus, the claim is established. 



min 



2n— l _ l _ l 

fc = 



n |ay-iw|-< 
3=1 

2n— 1 _1 _1 

E M 2 C k M 2, 
k = Q 

3 = 1 

/ M- 1 2 E C fe i W fe 
fe=Q 

M _1 2 J] 1 Ckiu' 

fc=0 



>-/3 



□ 



The search for minimum over all w £ R may in practice be confined to a small 
interval as 

/ 2n-l \ 

/ Af- 1 J2 i k C k u k \ 

k=0 



lim 

— > ±oo 



0. 



The above construction is illustrated by the following example. 
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Example 6.3. Consider the problem from Examples 12.11 and 12.21 of finding 
F(s) so that 



F(s) ( I) — > (75 

3/ V50 



From (|1.10[) we have that here 



'liiii 1 

1 2 4 8 16 32 
1 3 9 27 81 243 
1-11-11 -1 
1-2 4-8 16 -32 
.1 -3 9 -27 81 -243 




and thus 

Pi(s) = -3s 4 + 34s 2 - 13 
and from flSTj) *(s) is as in ([7T2|) . Hence, for all (3 e R, 

(6.4) F 4 (s) = P 4 (s) + /3*(s) = -3s 4 + 34s 2 - 13 + /3(1 - s 2 )(4 - s 2 )(9 - s 2 ), 
is an interpolating polynomial in Even. Next note that P 4 (s) is not in GPS, 

p 4 s )\s=^ = -( 3 ^ 4 + 34 ^ 2 + 13 ) 
namely in fact — P G C7P£. Next, from (|6.4j) one has that 

^4(s)| 8=to = (P 4 (s)+^(s)) u=to = (w 6 +14w 4 +49w 2 +36) + i(w 2 -l) 2 (w 2 +10). 
Thus, P 4 e C?P£ for all /3 > i. Furthermore, on the boundary, i.e. for (3 = i 



minRc (F 4 (s)| jt=iu) ) = P 4 (s)| s=± . = 0. 



Thus indeed one obtains all interpolating QV£ polynomials. 
Obviously, formally applying (|6.3|) leads to the same conclusion. 



Recall that in Example 12. 2[ through the simpler Lagrange approach we identified 
in (|2.4p only a subset of interpolating functions with /3 > | . □ 

7. A REFINEMENT 

We here address a case enabling us to refine the construction in the previous section 
by allowing in ()6.1[) M e P m (singular) and consequently obtain lower degree 
interpolating polynomials P(s). 

Recall that in subsection l5.1I Pfs) obtained, is an interpolating polynomial in Even. 
Thus on iR it is Hermitian. Assume that there exists a (constant) nonsingular 
matrix T so that 

(7.1) TP{s)T* = diag{P r (s), P m - r (s)} P m _ r E QVE, r 6 [0,m - 1]. 

(The case addressed in the previous section corresponds to r = m.) 

If r = 0, P(s) in (jl.lOp is already in GPS, it is a minimal degree interpolating poly- 
nomial and one can now employ (|6.2[) with arbitrary ^(s) from (|6.1[) and arbitrary 
/3 > 0. This is illustrated next. 

Example 7.1. Consider the scalar problem of finding all QPE polynomials P(s) 
so that 

*« (i)^(I 
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From fjl . 10[) we have that here 



■liiii 1 

1 2 4 8 16 32 
1 3 9 27 81 243 
1-11-11 -1 
1-2 4-8 16 -32 
,1 -3 9 -27 81 -243 





and thus, 



and from (|6.ip 

(7.2) *(*) = 

As we here have P(itu) = uj 2 
polynomials are given by 

(7.3) F(s) = P(s) + = 



P{s) = -s 2 + 5 
(1- S 2 )(4- S 2 )(9- S 2 ). 

f 5 for uj G M, P G QVZ and all interpolating 



+ 5 + /3(l-s 2 )(4-s 2 )(9-s 2 ) 



P > 0. 



It is of interest to mention that employing the simpler Lagrange approach described 
in Example 12.21 would identify the subset of interpolating polynomials F(s) in (|7.3|) 



with /3 > j^j. In particular, it would have failed to find the minimal degree inter- 
polating QVZ polynomial corresponding to ft = 0. □ 

Consider now the case where in (|7.ip 

r € [l,m- 1]. 

Then, conforming to P(s) in ()7.1|) . \&(s) in ()6.ip can be constructed, with the same 
T, so that 

T*(s)T* 



s)diag{M r , M m _ r }(x* + s) 



where M m _ r 6 



j=i 

_ r is arbitrary, including zero. 



M r G 



One can now proceed as before. Namely, fix M r G P r , where ||M r ||a = 1. Then /3 
in is given by 




iu>\ 



\i {M- l P r {iu)) 



Example 7.2. Consider the following two dimensional problem: Find all F{s) in 
QVZ so that 

F(l) = diag{-35, 9} 
_F(2) = diag{-20, 0} 
F{3) = diag{45, 25} 

From (ll.lOp we have that here 

/diag{-35 9} \ 
diag{-20 0} 
diag{45 25} 
diag{-35 9} 
. diag{-20 0} . 
\ diag{45 25}/ 



( J 2 


h 


h 


/2 h 


h \ 




/Co 


2/2 


4/2 


8/2 16/2 


32/ 2 1 




' Ci 


h 


3/2 


9/ 2 


27/2 81/2 


243/2 




C 2 


h 


-/2 


h 


-I2 I2 


-/2 




C 3 


\ 12 


-2/ 2 


4/2 


-8/2 16/2 


-32/2 j 




I c 4 


-3/ 2 


9/2 


-27/2 81/ 2 


-243/2 / 


\C 5 



Thus Ci, C3, C5 vanish and 
C = diag{-36 16} 



C 2 = diag{0 - 



Ca — F 
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Namely, 

P(s) = J 2 s 4 +diag{0, -8}s 2 +diag{-36, 16} 
= diag{s 4 -36, (4-s 2 ) 2 }, 

is a minimal degree interpolating polynomial in Even. Using (jl.3p the McMillan 
degree is 10. 

Note now that P(s) is of the form of (|7.1|) with T = I<z and r = 1. 
Now from (|6.1j) one can construct ^(s) with M £ P2. Indeed taking 

*(a) = (1 - s 2 )(4 - s 2 )(9 - s 2 ) (diag{/3, 0} + A) A £ P 2 , 

guarantees that 

F(s) = P(s) + *(s) 
is a <5"P£ interpolating polynomial for all 

P>P = 1. 

Moreover for A = one obtains 

F(s) = diag{s 4 -36+/3(l-s 2 )(4-s 2 )(9-s 2 ), (4-s 2 ) 2 } 

diag{-/3s 6 + (14£+l)s 4 -49,3s 2 +36(£-l), (4-s 2 ) 2 } 
= s 6 diag{-£ : 0} + s 4 diag{14,3+L l} + s 2 diag{-49,3, -8}+diag{36(/3-l), 16} 

where the McMillan degree is only 12 (see(fO])l, For /3 > 1 indeed F £ QVZ. □ 

8. A NON-CONVEX SET: {/-GENERALIZED POSITIVE EVEN POLYNOMIALS 

We here address the the interpolation problem where T is not convex. We 

modify the above "recipe" accordingly. As a test case, we take the set vQV£ 
described in (jl.9p . First, we show that this set is indeed not convex. 

As already mentioned, ±QV£ are subcics of Even, see e.g. [TTJ Section 5]. However, 
for v £ [1, in — 1] the set vQVZ in (|1.9j) is an invertible cone, but not convex. This 
can be illustrated even by constant 2x2 matrices. Indeed if one takes 

diag{-l, 4} = A = _R a diag{-l, l}i?* R a = diag{l, 2} 
diag{4, — 1} = B = _R. b diag{-l, l}R* b R b = (° g) 

then A + B = 3/2 is not a vQVZ function of the form of (|1.9[) . 

The fact that the set vQVE is not convex does not allow us to employ the Lagrange 
approach from Section [3J However, we can adapt the recipe from Section [TJ Here 
are the details. 

Given a reduced data set x\ , ... , x n £ C, the corresponding Y\ , ... , Y n £ C mxr ' 
and is, v £ [1, m— 1]. Find (all) F(s), low degree m x m-valued interpolating poly- 
nomials F : C -> C mxm , i.e. 

F(s) = G(s)diag{-/„ , 7„,-„}G*(s) t>e[l, m-1] G(s) polynomial 
F(^) = Y 3 i=l, ... , Ti. 

In order to have the problem feasible one needs to assume that the data satisfies 

Xj +xl=0 ==> Yj=Y k n>k>j>l 

Xj ei«. =>■ Y,-=T 3 diag{-/„ , 7 m _„}T* j=l, ... , n, 

for somc0 v £ [1, m — 1] and some Tj £ C mxm . 

''Recall, if 1/ = or v = £ we are essentially back to the QV€ case of subsection \E\ 
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Substituting in subsection 15.21 A = J = I m one obtains from Ijl.lOjl P(s), the 
minimal degree interpolating polynomial in Even. 

To construct the neutral polynomials ^f(s) substitute in (|4.ip M — i?diag{— /„ , 7 m _„}i?* 
with v € [1, 77i — 1] and above, to obtain 

n 

(8.1) := Y[( Xj - s) {Rdmg{-I v , I m - V }R*) (x* + s) Re C mxm , 

i=i 

with R parameter. 

By construction, with the above P(s) and ^(s), for all 8 G R 

(8.2) F(s)=P( S ) + /3tf(s), 

is an interpolating polynomial (jl.lj) in Even. However, as we already remarked, the 
set vQVE not convex. Thus one needs to justify the existence of 8 so that F(s) in 
is a vQVE interpolating polynomial for all 8 > 8. 



The idea relies on the following fact, here formulated in the framework of matrix 
theory. 

Lemma 8.1. Let A V ,A^ £ c mXm Hermitian matrices. Assume that A<s, has ir 
and m — tt eigenvalues in C+ and<C-, respectively (Ay is nonsingular) . Then, for 
all 

P > \\A P \\\\A^\\ 

the matrix (A p + 8A^) has ir eigenvalues in C+ and m — it eigenvalues in C_. 

This may be deduced in several way^, e.g. Weyl's Theorem [43l Theorem 4.3.1]. 
A detailed proof in our framework, is given in Proposition 18.21 below where we 
estimate 8. 

Proposition 8.2. Let F(s) be as in (|8.2[) with P(s) from (|1.10[) and ^f(s) as in 
(|8.ip where R is given and nonsingular. 
Then, in ([HUD F € vQVE for all f3 > (3 where 



Xj ~ iuj\ 2 



2n-l 



R-'CkiR*)- 1 ^ 



k=0 



(8.3) 8 = max IT ., , 

If x j — i^j with u>j € R is an interpolation node, a whole neighborhood of U)j is 
excluded from the above max. 

Proof Indeed, note that F(s) in (|8.2p is in vQVE if there exist a polynomial G(s) 
so that 

for almost all s£iR (P(s)+£*(s)) =G(s)diag{ -J„ , I m _ v }G*{s) 



for almost all sGiM £ C fc s fc +/3 f] -s)fldiag{-J„ , J m _„}i?*(a;*+s) =G(s)diag{- J„ , i",„_„}G # (s). 
\ fc=o j = i / 

If s = Xj e iR is an interpolation node, then 

2n-l 

E ^^U-, = lj = 2idiag{-7„ , I m -u}T* 

fc=0 



^In operator theory this is formulated as having A v + Aq, invertible whenever Ay is invertible 
and H^" 1 !!" 1 > ||A P ||, see e.g. [49] Theorem 10.20]. 
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and 

n 

J\{xj - s)Rdiag{-I„ , I m - v }R*{x* + s) = 0, 
i=i 

so the condition is satisfied in a neighborhood of Xj. Hence, assume hereafter that 
s £ iM. is out of a neighborhood of interpolation points. Thus, the above condition 
may be written as having a rational function G(s) so that 

for almost all sGiR -*=2; h^diag{-7„ . /,„_„} =G(s)diag{-/„ . 7 m _„}G#(s) 

Now, this in turn is implied by, 



for almost all s£i]R /3> 



E R- 1 C k {R*y 

k = 



n (.,-.)(.*+.] 

Substituting s = icu, u £ R, establishes the claim. □ 



Note that the value of (3 in (|8.3j) depends on the choices of R in (|8.1j) and of the 
norm in (|8.3|) . 

9. FUTURE RESEARCH 

As it is often the case, this study opens the door for future research problem. Wc 
here mention a sample of them. 

(i) Recall that in (11.8[) we mentioned that F £ QVE if and only if it admits a 
factorization of the form F(s) — G(s)G#(s). 

In Sections |3] through [7] we presented a "factorization free" recipe for ob- 
taining all minimal degree interpolating QVE polynomials. 
It is now of interest, for a given minimal degree interpolating QVE polyno- 
mial F(s), to explore properties of its (pseudo) spectral factors G(s), see 
e.g. [HI Section 5.2], [21] and Section 19.3]. 

(ii) An idea we have used throughout the work is as follows. If P(s) maps 
x\, ... , x n to Yi, ... Y n and ^(s) vanishes at Xi, ... , x n taking 

(9.1) F(s) = P(s) + P*(s) 

yields F(s) mapping X\, ... , x n to Y\, ... Y n , for all /3 £ C. 

Let now J 7 * , J-p be two families of functions where J 7 ^ is a (convex) subcone 
of Tp. Assuming P £ Tp and ^> £ J 7 *, for /3 > "sufficiently large" F(s) 
in (|9.1j) is an interpolating function within the family Ty. 



In [5] we adapt this idea to interpolation by scalar rational positive functions 
(as in the classical Nevanlinna-Pick Interpolation) but where the nodes are 
within C_. It is shown that there always exist interpolating functions of 
degree equal to the number of nodes. Moreover, an easy-to-compute recipe 
of constructing these functions is introduced. 

(iii) Adapt the recipe in Section[T]for interpolation by QV (not necessarily even) 
polynomials. The framework of (|9.1I) still holds with P{s) mapping Xj to 
Yj and ^(s) a minimal degree QV polynomial vanishing at Xj. However, 
there are two basic differences: 
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(a) There is no restriction on the structure of P(s). Namely, it is no longer 
necessary to double its degree. 

(b) ^f(s) is no longer of the form of (|4.1j) . For example, for two points 

xi,x 2 e C + 

*(s) = {x l - s)(x 2 - s)(6xl + (1 - 6)x* 2 + s)M 

where the scalar 6, 6 G [0, 1] and the matrix M , (M + M*) e P are 
parameters. 

(iv) Adapt the interpolation scheme of this work to cope with tangential struc- 
tured matrix-valued polynomial interpolation. Namely substitute (jl.ll) by 
HQ) and HQ). 

(v) One may be interested in interpolation by polynomials whose symmetry on 
the imaginary axis is of a group type, e.g. unitary, J-unitary, contraction 
or J-contraction, in the spirit of e.g. [S], [31] • It is of interest to solve 
the interpolation problem in where J 7 is a family of matrix valued 
polynomials with this symmetry. 

(vi) In Section [1] we pointed out that the set of QVE functions is a subcic 
(Convex Invertible Cone) of Even. If one focuses, as in this work, on the 
respective subsets of QVE and Even polynomials, invcrtibility is no longer 
relevant. Hence, out of the cic structure, it is only the Convex Cone part 
that can be used (as we indeed did). 

In (|2.3p we recalled the notion of reverse of a polynomial. It is then easy 
to verify that if a polynomial is in QVE, so is the corresponding reverse 
polynomial, (12.3[) . In fact, it turns out that the set of QVE polynomials is 
a sub-Convex Reversible Cone of Even polynomials. 

It is of interest to explore the Convex Reversible Cone structure of QVE 
polynomials and then to try employ it to interpolation. 
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